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ABSTRACT 

We present high-resolution, optical images (BVI + Ha) of the multiphase interstellar medium (ISM) 
in the thick disks of the edge-on spiral galaxies NGC 4013 and NGC 4302. Our images from the Hubble 
Space Telescope (HST), Large Binocular Telescope, and WIYN 3.5-m reveal an extensive population 
of filamentary dust absorption seen to z ~ 2 — 2.5 kpc. Many of these dusty thick disk structures have 
characteristics reminiscent of molecular clouds found in the Milky Way disk. Our Ha images show the 
extraplanar diffuse ionized gas (DIG) in these galaxies is dominated by a smooth, diffuse component. 
The strongly-filamentary morphologies of the dust absorption have no counterpart in the smoothly 
distributed Ha emission. We argue the thick disk DIG and dust-bearing filaments trace physically 
distinct phases of the thick disk ISM, the latter tracing a dense, warm or cold neutral medium. The 
dense, dusty matter in the thick disks of spiral galaxies is largely tracing matter ejected from the 
thin disk via energetic feedback from massive stars. The high densities of the gas may be a result of 
converging gas flows. This dense material fuels some thick disk star formation, as evidenced by the 
presence of thick disk H II regions. 

Subject headings: dust, extinction - galaxies: individual (NGC 4013, NGC 4302) - galaxies: ISM 
galaxies: spiral - galaxies: structure - ISM: clouds 



1. INTRODUCTION 

The feedback of energy and matter from stars to the 
interstellar medium (ISM) plays a central role in the evo- 
lution of galaxies. It impacts the physical/thermal state 
of the ISM, affecting the ability of the gas to form stars 
and affects the galactic-scale star formation by ejecting 
material from the thin disk or altogether from a galaxy 
(e.g., Dekel & Silk 1986). Most disk galaxies with suffi- 
cient star formation rates (or star formation rate surface 
densities), show extended "thick disks" of warm inter- 
stellar matter extending ~ 2 kpc away from the galaxy 
mid-plane (Rossa & Dettmar 2003a, 2003b). The pres- 
ence of a thick disk ISM supported against gravity is 
a direct result of energetic feedback processes (although 
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some additional material may be contributed through in- 
falling matter; see Sancisi et al. 2008). Thus, the study of 
the statistical properties (e.g., Rossa & Dettmar 2003a, 
2003b; Howk & Savage 1999) and detailed physics (e.g., 
Haffner et al. 2009) of the thick disk ISM can help us 
understand the role that feedback plays in disk galaxies. 

The energetic processes that expel gas from the thin 
disk also act on the solid phase of the ISM, interstellar 
dust grains. This has been utilized to study the physi- 
cal structure as well as the phase structure of the thick 
disk ISM (Howk & Savage 1997, 1999, 2000; Alton et al. 
2000, Thompson et al. 2004). While the energetic expul- 
sion of matter from the disk can impact the dust-to-gas 
ratio through dust destruction, the fact that prevalent 
extraplanar dust is seen in ~ 50% of normal spiral galax- 
ies (Rossa & Dettmar 2003a; Howk & Savage 1999) and 
even in galactic superwinds such as seen in M82 (Kaneda 
et al. 2010; Roussel et al. 2010; Engelbracht et al. 2006; 
Hoopes et al. 2005; Scarrott et al. 1991) demonstrates 
that dust can survive the ejection from the thin disk. 
In fact, the detection of thick disk dust via the extinc- 
tion it produces against background starlight in direct 
broadband images of edge-on galaxies is the most ob- 
servationally efficient method for detecting extraplanar 
interstellar matter (Howk & Savage 2000). 

The first detailed discussion of the implications of ex- 
traplanar dust was in the Howk & Savage (1997) (here- 
after HS97) study of the nearby edge-on spiral galaxy 
NGC 891, although the presence of high-z dust clouds 
in this galaxy were noted earlier (Keppel et al. 1991; 
Sandage 1961). HS97 showed that the thick disk of 
NGC 891 was threaded with a network of dusty clouds 
and filaments observed to heights up to z ~ 2 kpc. The 
most remarkable implication of this work was that the 
individual dust-bearing clouds have very large masses, 
often in excess of 10 5 M Q . Howk & Savage (1999) (here- 
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after HS99) subsequently surveyed a group of nine nearby 
edge-on spiral galaxies, demonstrating that high- z dust 
is found in many galaxies and is correlated with the pres- 
ence of extraplanar diffuse ionized gas (DIG) seen in Ha 
imaging. Subsequently, Rossa & Dettmar (2003a, 2003b) 
showed that ~ 40% of a much larger sample of spiral 
galaxies exhibit extraplanar dust and verified the cor- 
relation between the presence of extraplanar dust and 
DIG emission. HS99 and Howk & Savage (2000; here- 
after HS00) argued that the gas traced by extraplanar 
extinction features represented a much denser medium 
than that seen in Ha emission, with densities an order of 
magnitude or more higher than DIG material (see also 
Keppel et al. 1991). The extraplanar dust clouds and fil- 
aments, which are not spatially correlated with the DIG, 
are only visible because they have significantly higher ex- 
tinction - and likely higher column densities and particle 
densities - than their surroundings. HS00 hypothesized 
that the dust structures trace a cold, neutral medium 
(CNM) in the thick disk based on the large column den- 
sities, particle densities, and masses of the clouds. 

This conclusion was motivated in part by a detailed 
comparison of the extraplanar dust filaments and DIG 
emission by HS00, which demonstrated that they were 
tracing independent structures in the thick disk ISM of 
NGC 8910 Following the study of NGC 891, Thompson 
et al. (2004; hereafter THS04), used the Hubble Space 
Telescope (HST) to investigate the extraplanar dust in 
the edge-on galaxy, NGC 4217. These authors indicate 
that in addition to the smoother distribution of the Ha 
emission morphology compared to the extraplanar dust 
structures, they find little spatial correlation between 
areas of enhanced Ha emission around regions of con- 
centrated OB associations and high-z dust structures. 
Rossa et al. (2004) similarly used images from (HST) to 
show the ISM traced by extraplanar dust was distinct 
in morphology from the DIG in NGC 891, while Rossa 
& Dettmar (2003a, 2003b) argued for a similar arrange- 
ment in their larger sample of galaxies. The dusty ex- 
traplanar filaments are physically distinct from the DIG, 
even though the presence of these two tracers of the ex- 
traplanar ISM is strongly correlated (Rossa & Dettmar 
2003b; HS99). The small scale structures in the HS00 
and subsequent Rossa et al. (2004) imaging, the high 
column densities (Njj 10 21 in many cases), large cloud 
masses, and high particle densities (bh > 1 — 25 cm~ 3 ; 
HS00, Howk 2005) are all features in great contrast to the 
other observed phases of the thick disk ISM and point to 
a denser, cooler medium than the DIG. Thus, HS00 con- 
cluded that all of the thermal phases of the ISM found 
in the thin disks of spiral galaxies were also present in 
their thick disks. 

The only galaxy for which a detailed comparison of the 
Ha (DIG) and extraplanar dust morphologies has been 
undertaken is NGC 891. While the structures traced 
by these two probes of the extraplanar ISM were almost 
completely distinct in NGC 891, the universality of the 
result is not settled. In particular, Ha images of the 
DIG are often taken at relatively low spatial resolution, 
in part due to the faintness of the emitting gas, while 
the dust structures are lost to confusion at similar res- 

6 Keppel et al. (1991) reached the same conclusion on the basis 
of densities estimated for the DIG and extraplanar dust features. 



olution. NGC 891 has one of the brightest known Ha- 
emitting DIG layers (e.g., Rand 1996, Miller & Veilleux 
2003) and has more prominent filaments within the DIG 
than many galaxies (Heald et al. 2007). For galaxies 
with lower star formation rates, where the brightness of 
the DIG and the prominence of filaments diminishes, are 
there corresponding changes in the extraplanar dust? In 
this work we present high-resolution images of extrapla- 
nar dust and DIG in the edge-on galaxies NGC 4013 and 
NGC 4302. Our broadband images include observations 
from HST that elucidate the smallest scale structures in 
the dusty thick disk, while our Ha images are among the 
highest spatial resolution acquired from the ground and 
have good continuum subtraction. We use these images 
to extend the detailed comparison of extraplanar dust 
and DIG properties to these two galaxies, which have 
fainter DIG emission and a smaller fraction of the emis- 
sion arising in filaments (i.e., a larger proportion com- 
ing from the smooth DIG component). In the end we 
find no strong physical connection between the DIG and 
the structures traced by high-contrast dust absorption in 
these galaxies, similar to what is seen in NGC 891. The 
DIG and dust-laden clouds trace fundamentally different 
material in the thick disk. 

Our work is organized as follows. We discuss the ob- 
servations and data reduction in fJ5J In 331 we discuss the 
properties of the extraplanar dust and associated gas, 
while in $4] we describe the properties of the extraplanar 
Ha emission (DIG). In $5] we provide a direct compar- 
ison of the extraplanar dust morphologies and the DIG 
detected in Ha. We discuss in ^the origins, nature, and 
properties of the multiphase ISM in both NGC 4013 and 
NGC 4302 and summarize the major results of our work 
in §3 

2. OBSERVATIONS AND REDUCTIONS 

The work in this paper is based on optical imaging of 
two galaxies - NGC 4013 and NGC 4302 - whose proper- 
ties are listed in Table [TJ We use broadband and narrow- 
band imaging from several sources: 1) HST/Wide Field 
Planetary Camera 2 (WFPC2)) images of both galaxies 
in the (WFPC2) equivalents of the BVI bands to study 
extraplanar dust absorption at the highest resolution; 2) 
ground-based broadband images of NGC 4013 and NGC 
4302 taken with the WIYN 3.5-m and LBT 2x8.4-m tele- 
scopes, respectively, to study extraplanar dust at larger 
radial and vertical distances from the galaxy centers than 
allowed by HST. and 3) narrow-band Ha imaging of the 
extraplanar DIG in these galaxies, taken with the WIYN 
3.5-m, to study the connection (or lack thereof) between 
the ionized gas and material traced by dust absorption 
in the thick disks of these galaxies. We discuss the pro- 
cessing of these datasets in this section. 



2.1. HST/WFPC2 Observations and Processing 

We obtained broadband imaging observations of 
NGC 4013 and NGC 4302 with the Wide Field Planetary 
Camera 2 (WFPC2) aboard HST. WFPC2, described 
fully in Biretta et al. (2002), consists of four cameras. 
The three Wide Field Cameras (WFCs) together cover 
an "L" -shaped region 150" x 150" with intrinsic spatial 



Fig. 1. — Three color (BVI) composite mosaic of NGC 4013 created from the HST and WIYN imaging of this galaxy. The field of view 
of this image is ~ 5' X 5' on a side. North is up, east is to the left. Dusty thick disk clouds are seen as dark patches of extinction off of the 
galactic plane. The ffi>T/WFPC2 images of this galaxy cover the central and northeastern section of the galaxy but do not include much 
of the southwestern part of the galaxy. The bright emission in the plane just northeast of the bulge is a foreground star. There are slight 
edge artifacts due to the merging of the HST and WIYN images. 



sampling of O'.'l per pixel (f/12.9). The Planetary Cam- 
era (PC) covers a 34"x34" square field with 0'.'046 per 
pixel (f/28.3). 

Our observations, summarized in Table [TJ were ac- 
quired as part of GO program 8242 (PI: Savage). We 
imaged the majority of the galaxy NGC 4302 in one 
pointing, while we used two separate pointings to cover 
a significant fraction of NGC 4013. For NGC 4302 and 
one of the NGC 4013 pointings, we acquired eight ex- 
posures in the F450W filter and four exposures in each 
of the F555W and F814W filters. We obtained four 
images in each of the three filters for the secondary 
NGC 4013 pointing. The F450W, F555W, and F814W 
filters roughly correspond to the Johnson-Cousins B, V, 
and I bandfl 

7 We chose the F450W filter as the B-band equivalent because 
it has a higher over-all throughput than the F439W filter; the 
transformation to the standard Johnson-Cousins system is well- 



The images acquired for each pointing were split 
equally between two dithered sub-positions, with a linear 
dither step of 3.5 WF CCD pixels between the positions. 
We used the Dither II package and its drizzle algorithms 
(see Fruchter et al. 1998 and Fruchter & Hook 1997) to 
combine the dithered images onto a finer grid, following 
THS04. The resulting cosmic ray-cleaned images have 
an angular sampling of 0'/05 per pixel. 

Figures [1] and [5] show composite color images of 
NGC 4013 and NGC 4302, respectively, created from the 
WFPC2 BVI images. The WFPC2 images have been 
combined with ground-based observations in these im- 
ages. The ground-based images are more sensitive to 
faint, extended structures (i.e., the faint thick disk stel- 
lar emission) and probe a larger field of view. For NGC 
4013, the ground-based data are from the WIYN 3.5- 
m telescope and are discussed in H2.2I The ground-based 

behaved (see Holtzman et al. 1995). 
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Fig. 2. — Three color (BVI) composite mosaic image of NGC 4302 created from the HST and LBT imaging of this galaxy. The field of 
view of this image is ~ 8' X 5' on a side. North is to the right, east is up. Dusty thick disk clouds are seen as dark patches of extinction off 
of the galactic plane. The ffi>T/WFPC2 images of this galaxy cover the central region out to R ~ 6 kpc. NGC 4302 shows no signature 
of interaction with the neighboring galaxy, NGC 4298, seen in the bottom of the composite image. There are slight edge artifacts due to 
the merging of the HST and LBT images. 



observations of NGC 4302 were obtained with the 2x8.4- 
m Large Binocular Telescope (LBT) and are described in 



2.2. WIYN Observations and Processing 

To study the distribution of the diffuse ionized gas in 
these galaxies, we acquired narrow-band Ha and Gunn r 
images of NGC 4013 and NGC 4302 at the WIYN 3.5-m 
telescope at Kitt Peak National Observatory. These data 
were collected on 2004 April 18 and 19 (UT) under non- 
photometric conditions using the WIYN Tip-Tilt Mod- 
ule (WTTM), although without the tip-tilt functionality 
enabled. A log summarizing the WIYN observations is 
given in Table [TJ including the total exposure time in 
each filter and the seeing-limited resolution in the final 
images for eachgalaxy. 

The WTTM □ is an optical re-imaging system working 
at f/7.43. This system feeds a 2048 x 4096 EEV CCD 
with 13.5 /im pixels. Each pixel subtends 0'.'1125 on the 
sky giving a 3. '8 x 4f 7 field of view that is vignetted at 
the corners. The data acquisition, calibration, and re- 
duction of the images were done in a manner similar to 



8 The WTTM instument d escribed 
http : //www . noao . edu/wiyn/WTTM_manual . html 



here: 



that described in HS00. We note that because the field 
of view of the imager is relatively small compared with 
the galaxies themselves, the sky background subtraction 
is somewhat uncertain. None of our conclusions are sen- 
sitive to this problem; however, it precludes us from, e.g., 
properly deriving scale heights of the diffuse ionized gas 
in these galaxies. 

We derived the astrometric plate solution for our im- 
ages using a grid of stars whose coordinates were mea- 
sured in the Digitized Sky Survey for the NGC 4302 
observations. These agree well with the iKT-provided 
plate solution for the WFPC2 data. The rms error in 
using the plate solution to derive coordinates is always 
significantly less than l'.'O. For the NGC 4013 observa- 
tions, this treatment of the data left a systematic offset 
between the WIYN and HST data. Due to a lack of DSS- 
detected stars in the small WFPC2 field, we used stars 
and galaxies in common between the WFPC2 and WIYN 
images to bring the WIYN data onto the same reference 
frame as the HST data. The BVI imaging of NGC 4013 
was obtained with the previous generation WIYN Im- 
ager, the S2KB Camera. The properties and reductions 
of these data are as those discussed in HS00. These data 
are also summarized in Table [TJ 

The principal purpose of our ground-based WIYN ob- 
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servations is to study the distribution of Ha emission 
from the DIG in these galaxies. The WIYN W15 filter 
was used as an on-band Ha filter. This filter, centered at 
6569A with a FWHM of ~ 73A, contains emission from 
the Ha and nearby [N 110 lines plus a stellar continuum 
contribution. We used the Gunn r observations to es- 
timate the stellar continuum and subtract it from the 
narrow-band images, scaling the rband images so that 
faint Milky Way foreground stars were completely re- 
moved from the Ha images after continuum subtraction. 
Non-saturated foreground stars and background galaxies 
are well subtracted in our final line images. We found no 
evidence for a need to adopt a different scale factor for 
the bulge of the galaxies compared with the outer disks 
of the galaxies (e.g., Rand 1996). The Gunn r filter used 
for the continuum estimate encompasses the Ha line it- 
self and may lead to slight over-subtraction of the line 
emission. This has little impact on the current paper, in 
which we limit ourselves to morphological considerations. 

2.3. LBT Observations and Processing 

We obtained broadband UBVI images of NGC 4302 
with the Large Binocular Cameras (LBCs) on the 2x8.4- 
m LBT under non-photometric conditions on 2008 May 
03. A log summarizing the LBT observations is given in 
Table [T] including the total exposure time in each filter 
and the seeing-limited resolution in the final images for 
NGC 4302. The LBT uses two co-pointed 8.4-m mir- 
rors to feed a pair of nearly-identical LBCs. The LBCs 
are described in Giallongo et al. (2008). The two cameras 
operate simultaneously at prime focus, one optimized for 
blue wavelengths and the other for red (both can provide 
V-band imaging, the approximate cross-over in sensitiv- 
ity). Each camera provides a ~ 23' x 23' field of view 
using a four CCD mosaic with 0'.'23 pixels. Inter-chip 
spacings are filled in by dithering the telescope over sev- 
eral exposures. A log of our observations is given in Table 
[Q We used the LBC-Blue camera to image NGC 4302 in 
the U and B filters, while the LBC-Red camera provided 
coverage of in the V and I filters. We use twilight sky 
exposures to perform the flat field corrections. 

Reduction of the LBC images requires combining a se- 
ries of dithered images, taking into account the distor- 
tions that occur over the very wide field of view. We 
largely follow the techniques described in Sand et al. 
(2009) using modified versions of scripts from B. Werner 
et al. (in preparation). The data are processed in two 
steps. The first uses standard reductions from the IRAF 
package mscred to trim and remove the overscan re- 
gion, subtract bias frames, and apply flat fields (derived 
from twilight flats). Subsequently, we used the SCAMP 
(Bertin 2002) and SWARP (Bertin et al. 2002) soft- 
ware packaged to produce the final images. SCAMP 
was used to derive the astrometric solution for each of 
the chips in each image. The final solution is based on 
cross-matching sources identified in the LBT images us- 
ing SExtractor (Bertin & Arnouts 1996) with an astro- 
metric catalog derived from DR7 of the Sloan Digital 

9 While the images contain both Ha and [Nil] emission, for 
brevity's sake we will hereafter refer to our images as He? images. 
The reader should be aware that the bandpass also contains the 
[Nil] transitions. 

10 http://www.astromatic.net/ 



Sky Survey (SDSS-DR7; Abazajian et al. 2009). This is 
a necessary step due to the distortions over the full field 
given the //1.14 focal ratio. The absolute astrometric 
calibration has an rms dispersion compared with SDSS 
of ~ 0'.'2 across the entire field of view. The individual 
images were subsequently resampled and coadded using 
SWARP. 

3. THE DUSTY THICK DISKS OF NGC 4013 AND 
NGC 4302 

Visual inspection of Figures Q] and [2] reveals the pres- 
ence of strong, patchy extinction in the thick disks of 
these two galaxies. Filaments and clouds of absorbing 
material are seen to nearly the radial edge of the galax- 
ies' optical disks and to z-heights approaching the largest 
distances to which stellar light can be seen in the rel- 
atively shallow WFPC2 images, although the ground- 
based LBT and WIYN images reveal the dust extinction 
extends to lesser heights than the stellar thick disk. This 
extraplanar extinction is a direct tracer of the thick disk 
ISM in these galaxies (e.g., Howk 2005). In order to be 
seen, these dusty regions must have higher column densi- 
ties than their surroundings with these filaments having 
excess extinction compared with nearby sight lines, pro- 
viding significant contrast between these clouds and their 
local environment in our images. Smoothly distributed 
dust does not produce such contrast; thus a diffusely dis- 
tributed component of dust may exist undetected in the 
thick disks of these galaxies. That the column densities 
of the observed dust-bearing clouds are higher than their 
surroundings almost certainly implies the dust features 
seen in our images have significantly higher volume den- 
sities than their surroundings, as well. 

Figure [3] shows the fiSTand WIYN V-band images of 
NGC 4013. The top panel shows the HST V-band data 
and the middle panel shows the unsharp masked version 
of the top panel. The display of the galaxy in the top 
panel is not able to show the absorbing clouds against the 
bulge and further out in a single display, making the un- 
sharp masking necessary. The unsharp masked versions 
of our images remove the large-scale gradients in the un- 
derlying stellar light, revealing structure on scales smaller 
than the smoothing kernel (FWHM < 350 pc) (see HS00 
and THS04 for detailed discussions of this approach). We 
only use these images for display and have taken care 
to assess and minimize masking-related artifacts. Both 
HST images include insets showing the central bulge re- 
gion with a stretch that shows the dust absorption more 
clearly. The bottom panel shows the unsharp masked 
display of our final WIYN data. The WIYN image gives 
almost full radial coverage of the galaxy, revealing extra- 
planar dust structures extending along nearly the entire 
radial length of the disk of NGC 4013. Extraplanar dust 
clouds can be identified against the background stellar 
light z ~ 2.0 kpc from the midplane. Those dust clouds 
at the largest heights show less prominent substructure 
and generally have lower ay values, but they are still well 
detected in our WIYN images. We see no evidence for 
disparities in the amount of extraplanar dust comparing 
opposite sides of the plane. 

Similarly, Figure H shows the HST and LBT V-band 
images for NGC 4302. The top panel shows the HST 
V-band data and the middle panel shows the unsharp 
masked version of the top panel. The HST images cover 
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Fig. 3. — HST & WIYN V-band images of NGC 4013. The top panel shows the HST V-band image; the middle panel shows the unsharp 
masked version of this image. The insets in the top two panels are stretched to show the bulge region in greater detail. The area of the 
insets is indicated by dashed lines. The bottom panel shows the unsharp masked version of the WIYN V-band image. The dashed box 
in the WIYN V-band indicates the extent of the HST image shown in the upper panels. All image displays are inverted so that areas of 
dust extinction appear lighter than their surroundings. The HST images show dense, narrow, filamentary-like extraplanar dust structures 
at high resolution. The WlYN image show the extensive dense dust structures and complexes extending to heights z ~ 2 kpc beyond the 
midplane. 



~ 6 kpc projected radial distance from the center on ei- kpc from the midplane, detected in the LBT images, 

ther side, while the LBT images show ~ 13 kpc radial The FHWM of the Gaussian used to produce the un- 

span of the optical disk and more. The dashed box in sharp masked images for NGC 4302 were the same as in 

the LBT V-band indicates the extend of the HST image NGC 4013. 
shown in the upper panels. Extraplanar dust clouds can 
be identified against the background stellar light z ~ 2.5 




Fig. 4. — HST & LBT V-band images of NGC 4302. The top panel shows the HST V-band image; the middle panel shows the unsharp 
masked version of this image. The bottom panel shows the unsharp masked version of the LBT V-band image. The dashed box in the LBT 
V-band indicates the extend of the HST image shown in the upper panels. All image displays are inverted so that areas of dust extinction 
appear lighter than their surroundings. The HST images show a large number of extraplanar dust structures at high resolution. The LBT 
image shows extraplanar dust-bearing clouds in the thick disk along nearly the entire radial length of the galaxy and is more sensitive to 
small extinction clouds. 



3.1. Quantifying Dust Extinction 

Using the HST V-band data for our measurements (top 
panels of Figure [3] and El and following HSOO and THS04, 
we provide spatially-averaged estimates of the apparent 
opacities of a small number of extraplanar dust struc- 
tures. From these we derive the "apparent" extinction, 
column densities, particle density, and masses of these 
dust-bearing structures. We define the "apparent" ex- 



tinction, ax, for a given wave band A as 

a x = -2.5log(S dc ,x/S bg ,x); (1) 

here Sd c ,\ is the average surface brightness measured to- 
ward a dust cloud (dc), and Sb g ,\ is the surface brightness 
of the local background (bg) . Defined in this way, the ap- 
parent extinction is a lower limit to the true extinction, 
Ax, because Sdc.x includes both extincted starlight from 
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behind as well as unextincted light emitted in front of the 
feature. One can, in theory, use a\ for multiple wave- 
lengths with an assumed extinction curve to estimate 
the fraction of light originating in front of the dust fea- 
ture (e.g., Gallagher & Hunter 1981, Knapen et al. 1991, 
HS97). In this case, the fraction of Sdc,x arising in front 
of the clouds is x and the true extinction is a function 
of ay and x. However, HS00 found that this approach 
does not always lead to self-consistent results, likely due 
to the effects of spatial averaging of probing differing 
depths into clouds at different wavebands and of scatter- 
ing. HS97 (who used a narrower range of wavelengths) 
provided evidence for a significant amount of stellar light 
arising in front of the dust clouds in NGC 891. 

Here we follow HS00 and THS04 in assuming Ay > ay 
(equivalent to assuming x > 0). This has the advantage 
that our calculations of the physical properties of the 
structures seen in our images will provide lower limits to 
the true properties. HS00 discuss the degree to which 
Ay departs from this assumption for several values of 
x. For example, for x = 0.25, the V-band extinction 
Ay w 1.5 a v . We do not report the extinction in the 
other bands. 

For each dust feature we measure the ratio of sur- 
face brightnesses Sd c ,xl Sbg,x using intensity distributions 
measured vertically or horizontally, with respect to the 
disk, averaging the results for a given cloud. We es- 
timated the background surface brightness appropriate 
for each feature by fitting a spline to the intensity dis- 
tribution of the star light in each cut. This technique 
is different than that used in our previous works where 
we were more concerned with the impact of the smear- 
ing due to the seeing. The current measurements will 
average the ay values over more of the clouds than those 
earlier measurements, leading to somewhat smaller ex- 
tinction values. 

3.2. Physical Properties of Individual Dust Structures 

Table Q] summarizes the physical properties of several 
cloud complexes in the thick disks of the two target 
galaxies from the HST V-band data. Figures [5] and |6] 
identify the clouds discussed in Table [TJ they are dis- 
played in the unsharp-masked WFPC2 images to better 
show the small-scale features of the selected clouds. Sev- 
eral measured quantities are given, including their posi- 
tions, physical dimensions, height above the midplane, 
and ay values. Estimates of the column densities, den- 
sities, and masses associated with the clouds are also 
given, as outlined below. In this work we have favored 
small, relatively isolated structures in order to study the 
densities implied by the small-scale structure observed in 
the WFPC2 imaging. The choice to work with isolated 
clouds tends to favor lower extinction clouds at relatively 
large z. The smaller sizes (to emphasize high densities) 
tend to favor lower masses. 

The values of ay in Table [1] are in the range 0.1-0.3 
mag. These values represent spatial regions of the clouds 
lower than the values appropriate for the cores of the 
structures. Given the values of ay are also lower than 
the true extinction, realistic values of ay for much of the 
clouds are likely significantly larger than those values 
given in Table [I] 

Transforming the apparent extinctions discussed above 
into physical properties of the clouds requires an as- 



sumption about the gas-to-dust relationship in the thick 
disk. Following earlier works, we assume the clouds fol- 
low dust to gas relationships appropriate for the disk 
of the Milky Way (HS97, HS99, HS00). This is a nec- 
essary step for estimating the physical conditions of 
the dust features from the measured apparent extinc- 
tions, one whose consequences are discussed at length in 
HS00. We adopt the average relationship between to- 
tal hydrogen and the color excess E(B — V) of diffuse 
clouds in the Galactic disk as estimated by Bohlin et 
al. (1978): JV(H I + R 2 )/E(B -V) = 5.8 x 10 21 atoms 
cm~ 2 mag -1 . To relate this to the extinction, we assume 
Ry = Ay/E(B — V)~ 3.1, roughly the mean for diffuse 
clouds in the Galactic disk (e.g., Cardelli et al. 1989). 
This gives the total hydrogen (H I + H2) column density 
from the apparent extinction: 

N H (cm- 2 ) > 1.9 x 10 2 V. (2) 

The inequality arises since ay < Ay. Since both 
NGC 4013 and NGC 4302 are relatively massive spiral 
galaxies, the assumption of a Milky Way gas-to-dust ratio 
will probably yield a reliable estimate of the associated 
gas column density. However, if significant dust destruc- 
tion has occurred in any particular feature, the derived 
Nh will be an underestimate, while if the dust has been 
separated from the gas then Nh will be an overestimate. 
We favor scenarios in which little dust destruction occurs 
in moving the dust from the thin to thick disk or where 
some destruction has occurred in which the dust is sepa- 
rated from the gas. We thus quote the column densities 
as lower limits. 

We give lower limits to the masses of for the features 
in Table [TJ calculated as M ~ fimnh t h Adc, where Ad c 
is the projected area of the dusty clouds. We assume 
fi = 1.37 to correct for the contribution to the mass from 
helium and the heavy elements. The definition of the 
cloud areas is a critical and highly-subjective component 
of the mass estimates. In the present work we have fa- 
vored small clouds and have been conservative in our 
derivation of the projected areas. The implied masses 
should be treated as lower limits. 

In addition to the columns and masses, we also in- 
clude a very crude density estimate in Table [1] This esti- 
mate assumes the clouds we see are largely cylindrically- 
symmetric such that nu ~ Nh/Ax, where Ax is the 
width of the minor axis of the clouds. The only support 
for this assumption is the general prevalence of filamen- 
tary absorbing structures seen in our images, although 
this impression could be affected by the general confu- 
sion caused by overlapping clouds along the line of sight. 
In general the morphologies of the absorbing clouds are 
complex, but we move forward using this simplistic as- 
sumption for estimating the densities. 

Keeping in mind the uncertainties discussed above, 
the properties of the clouds summarized in Table [1] are 
nonetheless intruiging. The features all have estimated 
columns Nh > 2 x 10 20 cm -2 , simply as a matter of se- 
lection (it is difficult to identify clouds with a v < 0.1). 
We infer cloud masses > 10 4 M Q for all but the smallest 
two clouds in Table [TJ and the densities implied for those 
clouds included in the table are of order tih ~ 1—4 cm -3 . 
Additionally, we examined some of the darkest narrow 
filaments in the WFPC2 images located in regions too 
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Fig. 5. — HST unsharp masked V-band image of NGC 4013 with select dense extraplanar clouds outlined and labeled. The properties 
of these clouds are summarized in Table [T] These particular absorbing structures were chosen based on their location, apparent opacity, 
and size in an effort to emphasize the physical properties of smaller, dense structures that the high-resolution HST images reveal. These 
extraplanar clouds are located at heights z ~ 430 — 930 pc, have estimated column densities and masses TVh > 2 x 10 20 cm -2 , and M > 10 4 
Mq. Substructure in these clouds is seen to the limit of our resolution. The orientation of this image is as in Figure [3] Figure \E\ gives a 
close-up image of the gray outlined clump to the east below the bulge, cloud 115834.4 + 435651.9. 




Fig. 6. — HST unsharp masked V-band image of NGC 4302 with select dense extraplanar clouds outlined. The orientation of this image is 
as in Figure [2] The properties of these clouds are summarized in Table [T] These particular absorbing structures were chosen based on their 
location, apparent opacity, and size in an effort to emphasize the physical properties of smaller, dense structures that the high-resolution 
HST images reveal. These extraplanar clouds, located at heights z > 590 pc, have estimated column densities and masses Nn > 2xl0 20 
cm -2 , and M >10 4 Mq, respectively. Substructure in these clouds is seen to the limit of our resolution. 



confused to show well in Figures [5] and |6l so they are 
not included in Table [TJ Several of these clouds have 
column densities approaching or exceeding 10 21 cm -2 , 
giving mean densities nu ^8 — 10 cm -3 , which is ex- 
traordinary for thick disk clouds. The column densi- 
ties and densities measured for the thick disk clouds in 
NGC 4013 and NGC 4302 are similar to those seen in 
NGC 891 and NGC 4217 (HS97, HS00, THS04, Howk 
2005). The derived densities likely depend on the reso- 
lution with which one observes the structures, as higher 
resolution images not only smear less background light 
into the structure but also allow the identification of 



smaller-scale filaments. The masses given in Table Q] are 
smaller than have been reported in earlier works (HS99 
HS00, THS04), but that is largely a function of our newly 
adopted measurement technique and the choices we have 
made for emphasizing the smaller subcomponents rather 
than the largest complexes. Complexes of absorption 
with implied total masses well in excess of lO 5 Ajf0 cer- 
tainly exist in these galaxies. Overall, we do not see any 
clear evidence that the dust-bearing thick disk clouds in 
these galaxies are significantly different in character than 
those seen in NGC 891 (HS00) or NGC 4217 (THS04). 
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Fig. 7. — (a) A section of the bulge of NGC 4013 from the HST unsharp masked V-band images showing the large dust complex discussed 
in the text, (b) Close-up view of the large complex, which extends ~ 900 pc by 200 pc. This complex has a mean extinction ay = 0.04 
and total mass > 10 s Mq. (c) An enlarged view of two of the cometary-like clouds associated with the large complex. These are clouds 
115832.2 + 435644.6 and 115832.0 + 435644.1 in Table[l] Cloud 115832.2 + 435644.6 has a notable core-halo type structure, as emphasized 
in the relative contours shown on this display. The core of the cloud has ay = 0.36 and implied density njj > 7 cm -3 . The orientation of 
these images is as in Figure [3] 



3.3. Apparent morphologies of dust-bearing gas clouds 

Our images allow us to identify large-scale dust-bearing 
clouds in the thick disks of these galaxies and study their 
internal structure. Our high angular resolution HST 
images reveal substructure within the cloud complexes, 
down to the limit of our resolution, but only for the in- 
ner sections of both galaxies. Our LBT and WIYN data 
allow us to identify dust clouds at larger radial and verti- 
cal distances due to their greater field of view and depth 
compared with the HST images. Similar to clouds in the 
Milky Way, in which structure is seen on many scales, the 
thick disks of NGC 4013 and NGC 4302 contain narrow 
filamentary structures, large complexes of clouds (specif- 
ically in NGC 4302), as well as small, very dense clouds 
seen to our resolution limit. 

3.3.1. NGC 4013 

The extraplanar dust-bearing clouds in NGC 4013 are 
generally highly structured filaments, although we see a 
broad range of morphologies in both the HST and WIYN 
images in Figure [31 The typical dust filaments are ~ 100 
pc by a few x 100 pc in size (see Table [T] and Figure 
[3]). We note that the dust filaments are unlikely to be 
associated with a flare in the outer galaxy (see detailed 
discussion in HS99, HS00) or as a result of the warp in 
NGC 4013, which is perpendicular to our line of sight 
(Bottema 1996) and therefore would not give rise to the 
symmetrical distribution of clouds seen in our images. 

There is a strong preference for the highest- z absorbing 
structures (z > 2 kpc) in NGC 4013 to appear bent to- 
ward the northeast, especially on the eastern side of the 
galaxy (see bottom panel in Figure ■ Such an arrange- 
ment could occur if there is a rotational lag that increases 
with height (e.g. Heald et al. 2007). This is consistent 
with the sense of rotation of NGC 4013, with the south- 
west side receding from us (Bottema 1996). NGC 4013 is 
a member of the spiral-rich Ursa Major group of galax- 



ies. Ursa Major is considered a loose group with no cen- 
tral concentration and no sign of strong H I deficiencies, 
suggesting this group lacks a strong intracluster medium 
(Angiras et al. 2007). As a member of the Ursa Major 
group, NGC 4013 is relatively isolated. Thus the bending 
of filaments is not likely the result of movement through 
an intracluster medium in this case (Martinez-Delgado 
et al. 2009). We find no evidence for obvious supershell 
structures in NGC 4013, although identifying such struc- 
tures within the complex distribution of clouds might be 
difficult. 

We see a set of very large-scale complexes of dust 
above and below the bulge of the galaxy in Figure |31 
These large dust complexes, traced out to z ~ 1.8 kpc 
on either side of the bulge, show significant small-scale 
structure within them, including dense clouds similar to 
those cometary shaped clouds identified in HS00. Fig- 
ure [7] shows one of these complexes to the southeast of 
the bulge. This large complex (~ 900 pc x 200 pc) ex- 
tends from within ~ 500 pc of the plane, bending to- 
ward the east with increasing height above the plane. 
It shows several cometary-like substructures at heights 
z ~ 600 — 850 pc. Figure [T^b) provides a close-up image 
of this large complex, which has a mean extinction value 
of ay = 0.04 and total mass of order > 10 5 M Q . Figure 
[7(c) shows two of the the cometary-like clouds, clouds 
115832.2 + 435644.6 and 115832.0 + 435644.1 from Ta- 
ble [TJ Contours are included to guide the eye (but are 
based on the unsharp masked image and are not quanti- 
tatively useful). Cloud 115832.2 + 435644.6 has a head- 
tail type structure; the head of the cloud has ay = 0.36 
and implied density nn > 7 cm -3 . The association of 
these cometary structures with the larger complex sug- 
gests the complex is interacting with an ambient medium 
as it falls toward the disk through a corona or outflowing 
wind from the disk. 

All of the head-tail clouds we see in NGC 4013 seem to 
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be associated with larger complexes, like that shown in 
Figure [71 The head-tail clouds do not seem to be com- 
mon, although the clouds may still be interacting with an 
ambient medium. Such morphologies are reminiscent of 
those seen in some high- velocity clouds (HVCs) passing 
through the diffuse halo of the Milky Way. Putman et al. 
(2011) investigated the H I 21 cm emission line morpholo- 
gies of isolated HVCs, finding only ~ 35% of them show 
head-tail morphology, suggesting such morphologies are 
short lived. If the extraplanar clouds in our images are 
analogous to the HVCs, those with cometary morpholo- 
gies may be moving at larger velocities and interacting 
more strongly with the thick disk ISM. While the HVCs 
of the Milky Way do not appear to contain dust with 
the exception of the Magellanic Stream (Lu et al. 1998, 
Sembach et al. 2001, Richter et al. 2001, Fox et al. 2010), 
and their origins therefore are different from the thick 
disk dusty clouds studied here, the physics that dictates 
the head-tail morphology is likely quite similar. That is, 
the magnetohydrodynamical processes that shape those 
cometary structures likely do not strongly depend on the 
dust content of the clouds. 

Our HST images reveal a vast array of small-scale ab- 
sorbing features with dimensions approaching our resolu- 
tion limits. The smallest-scale dust structures are smaller 
than a few x 10 pc in size. We have measured the sizes 
and opacities of a number of these very small structures 
and estimate they have ay ~ 0.2 — 0.4 mag and njj > 8 
cm~ 3 . Indeed, each complex exhibits significant sub- 
structure. To demonstrate this we show in Figure [5] an 
expanded view of NGC 4013 cloud 115834.4 + 435651.9, 
located at z ~ 930 pc from the midplane. Prominent 
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Fig. 8.— A close-up view of cloud 115834.4 + 435651.9 in 
NGC 4013. The orientation of this image is as in Figure [3] This 
large cloud is located z ~ 930 pc from the midplane. The contours 
show ay = 0.2 to 0.6 in steps of 0.1 mag. As a single cloud we 
estimate the mass to be M >7 X 10 4 M®, although this depends 
on the definition of the cloud boundaries. This cloud shows promi- 
nent substructure. Two particularly dense dense subclouds at the 
top and bottom of the cloud have ay > 0.5 mag, implying column 
densities of Njj > 10 21 cm -2 and particle densities of rijj > 7 
cm -3 . 

substructures are clearly seen. The contours in Fig- 
ure [8] range from ay ~ 0.2 — 0.6. For comparison our 



spatially-averaged value for this cloud is ay ~ 0.22 (Ta- 
ble [J) . This demonstrates the impact of our averaging 
technique when making our measurements. Two partic- 
ular subclouds which are similar in size and opacity, have 
ay > 0.5 mag, implying column densities of Nh > 10 21 
cm and particle densities of nn > 7 cm~ 3 . These 
are lower limits for the column densities and imply that 
these dust structures have properties similar to Milky 
Way diffuse clouds, yet are located ~ 1 kpc away from 
the midplane. 

While the dust in highly-structured clouds is most 
readily apparent in our images, one can find very low- 
extinction, more diffuse material in the thick disk in our 
images (particularly the WIYN and LBT images, which 
have greater sensitivity to low surface brightness fluctua- 
tions). The lowest extinction clouds visible in our images 
of NGC 4013 have ay ~ 0.03 mag. If we couple that with 
a rough measure of the area occupied by these low col- 
umn density clouds (with very low contrast), the total 
mass implied in this galaxy is of order ~ 10 s M Q . This 
is likely similar to the total mass in the population of 
denser, more readily apparent complexes discussed above 
(see HS00). 

3.3.2. NGC 4302 

The extraplanar dust-bearing clouds in NGC 4302 are 
filamentary with some large dust complexes also notable 
in the thick disk in both the HST a,nd LBT images in Fig- 
ured! The LBT unsharp masked display (bottom panel 
of Figure |4|) shows extensive extraplanar dust-bearing 
clouds along the entire length of the plane detectable 
to z ~ 2.5 kpc. While large complexes of extraplanar 
dust are seen on both sides of the planes, the intensity 
of starlight as a function of z shows a clear asymmetry 
across the plane in the bulge region in NGC 4302. The 
apparent deficit of starlight at z ~ 1 kpc on the west 
side of the plane corresponds to a large diffuse complex 
of dust on this side of the disk. This diffuse complex is 
not easily seen in the HST unsharp masked images given 
its size and relatively little substructure, but is visible in 
the unsharp masked LBT image (bottom panel of Figure 
[4]). This extinction complex on the western side of the 
bulge area may be the cause for a measured asymmetry 
in Ha emission image as well (see 

A common characteristic in the majority of extrapla- 
nar clouds in NGC 4302 is the apparent bending of the 
large complexes toward the north at large z-height, out 
to z ~ 2.3 kpc. This may be connected with the motion 
of NGC 4302 through the Virgo Cluster, as discussed by 
Chung et al. (2007). These authors postulate the exten- 
sive tail of the H I disk extending along the northern edge 
of NGC 4302 is due to the galaxy's motion toward the 
south as it falls into the Virgo Cluster. Our images may 
be showing a physical manifestation of this movement 
through the intracluster medium in the shaping of the 
largest z dust features. We note that NGC 4298, pro- 
jected close on the sky to NGC 4302, is also located in 
the Virgo Cluster. NGC 4302 shows no sign of a strong 
interaction with NGC 4298. 

The ratio of thick-disk to thin-disk stellar scale heights 
is smaller in NGC 4302 than many other edge-on galaxies 
(as it is in the Ha; Rand 1996). Combined with the line- 
of-sight overlap of structures, this makes it difficult to 
cleanly identify large scale extraplanar dust structures. 
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Our images do not provide strong evidence for any ob- 
vious supershells or galactic chimney features extending 
from the midplane of NGC 4302. Figure [6] shows the un- 
sharp masked HST images of this galaxy with the clouds 
from Table[T]marked. As in NGC 4013, we see small-scale 
structures in the dusty clouds in NGC 4302 to the resolu- 
tion limit of our HST images (~ 10 — 20 pc). The implied 
densities of the small-scale clouds range from rtn > 2 to 
14 cm -3 for the narrowest, darkest thick disk clouds. 

Similar to the case of NGC 4013, there are diffuse, 
lower column density structures seen over a significant 
amount of the thick disk region of NGC 4302. These tend 
to have minimum extinctions ofay ~ 0.02 mag, implying 
a total thick disk mass in this more diffuse matter of 
~ 10 8 Mq. 

4. IONIZED GAS IN THE THICK DISKS OF 
NGC 4013 AND NGC 4302 

Extraplanar diffuse ionized gas (DIG) has been previ- 
ously studied in both NGC 4013 (Rand 1996, Miller & 
Veilleux 2003) and NGC 4302 (Pildis et al. 1994, Rand 
1996, Rossa & Dettmar 2000, Collins & Rand 2001, 
Heald et al. 2007). We have revisited these galaxies in an 
effort to achieve higher resolution, greater depth, and/or 
more consistent continuum subtraction a previous works. 
Our overall goal is to compare the Ha and dust mor- 
phologies to study their connection in these galaxies, 
which requires high resolution imagery. We detect thick 
disk Ha emission, with 0.8" resolution, in both galaxies. 

4.0.3. NGC 4013 

Extraplanar DIG at large ^-heights can be seen over 
nearly the whole of the star-forming disk (to R ~ 9 
kpc) in NGC 4013 (see the bottom-right panel of Fig- 
ure [9j which emphasizes the faintest Ha emission) . The 
extraplanar Ha emission in NGC 4013 is considerably 
fainter and shows much less prominent filamentary emis- 
sion than in NGC 891 (Rand 1996). Rand (1996) gives 
Ha emission scale heights of h z ~ 0.7 — 0.9 kpc (for the 
eastern side); however, he also notes significant contin- 
uum subtraction difficulties affect his images. Miller & 
Veilleux (2003) fit a scale height of h z ~ 1.5 kpc (for 
the "southeastern" side). Figure O shows our Ha im- 
ages of NGC 4013 in four panels; the top-left panel is 
an unsharp-masked WIYN V-band image to show the 
dust absorption at the same resolution as the Ha emis- 
sion while the bottom-left panel is an unsharp-masked 
Ha image to show the small scale features (as before, 
the smoothing kernel has FWHM< 350). Both right 
panels in Figure [9] are Ha emission, stretched to em- 
phasize the brightest emission (top-right panel) and the 
faintest emission (bottom-right panel). Our images are 
consistent with the description of the DIG morphology 
for NGC 4013 given by Rand (1996): diffuse with few 
filaments. Over the majority of the disk, the DIG mor- 
phology is smooth and lacks bright, narrow filamentary 
structures like those seen in the Ha emission of NGC 891 
and NGC 5775 (Rand 1996, HS00). What filaments exist 
in Figure [S] are relatively subtle, the brightest of which 
have peak intensities of ~ 10 — 20% of the local smooth 
DIG. Even the brightest filaments in NGC 4013 are com- 
paratively fainter than those in NGC 891, which have 
30% - 100% of the brightness of the local diffuse DIG 



emission (HS00). Given the faintness of the smooth DIG 
layer in NGC 4013, the filaments have a much lower ab- 
solute intensity as well. 

One of the most notable morphological features of the 
DIG in NGC 4013 is the large "H" -shaped structure em- 
anating from both sides of the bulge region (see Rand 
1996). This "H" structure consists of four large-scale fil- 
aments present at similar positions on either side of the 
bulge, at projected radial distances R ~ 200 — 300 pc and 
visible in our images to z ~ 2.2 kpc. Rand (1996) argued 
this feature was a relic nuclear (bipolar) outflow from 
an earlier nuclear starburstor AGN. Garcfa-Burillo et al. 
(1999) mapped the "H" structure in CO emission (see 
their Figure 3) and argued for a close association between 
the extraplanar DIG filaments and the extraplanar CO 
emission. Schinnerer et al. (2004) also mapped CO emis- 
sion from the plane of NGC 4013, connecting the "H" 
structure in a high-density CO ring in the disk. These 
works both invoke a superwind outflow as the origin of 
this structure, in agreement with Rand's conclusion. 

4.0.4. NGC 4302 

Figure [10] shows the continuum-subtracted Ha emis- 
sion in NGC 4302. The images of this galaxy are not 
as deep as for NGC 4013, but they have better resolu- 
tion than most previous Ha images of this galaxy (Pildis 
1994, Rand 1996, Rossa & Dettmar 2000, 2003a, 2003b). 
NGC 4302 has one of the faintest detectable DIG lay- 
ers (Rand 1996, Rossa & Dettmar 2000) Collins & Rand 
(2001) measure the emission scale height h z ~ 0.55 — 0.7 
kpc. Rand (1996) characterizes the extraplanar DIG 
morphology of NGC 4302 as "faint and diffuse" ; our im- 
ages are consistent with both of those characterizations. 
In our images the thick disk DIG is smoothly distributed 
within the first z ~ 2 kpc of the plane, but only within 
R ~ 7 kpc and 8.5 kpc radially from the bulge in north- 
ern and southern halves of the galaxy, respectively. The 
Ha emission abruptly decreases to undetectable levels 
at larger radii. We detect slightly brighter extraplanar 
DIG on the eastern side compared to the western side of 
the galaxy, most notably around the bulge region. This 
asymmetry is also seen in continuum star light and is due 
to enhanced extinction from large extraplanar dust com- 
plexes on the western side of the plane at z ~ 1 kpc (dis- 
cussed in §3.3). We find no bright plumes or filamentary 
structures. There are some very weak DIG filaments, the 
brightest of which measure only ~ 20 — 25% above the 
local (very faint) diffuse background. Thus, the DIG in 
NGC 4302 shows little to no structure or filaments, and 
what substructure is detected seems to be the result of 
dust extinction affects. 

Dettmar (1995) and Heald et al. (2007) found that 
the presence of Ha emission in NGC 4302 was weak 
and diffuse, especially in comparison to NGC 891, al- 
though both galaxies have prominent extraplanar dust. 
Our Ha images are the highest resolution available for 
both NGC 4013 and NGC 4302 and show both have 
smoothly distributed DIG layers with much less promi- 
nent filaments than the DIG of NGC 891 and even less 
than NGC 5775 (Rand 1996, HS00). Both NGC 4013 
and NGC 4302 also have comparably weaker Ha emis- 
sion than NGC 891 and NGC 5775, which are known for 
their bright Ha emission and have larger DIG emission 
scale heights (Collins & Rand 2001). Thus the galaxies 
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Fig. 9. — Comparison of WIYN V-band and continuum-subtracted Ha emission-line images of NGC 4013. The left panels shows the 
unsharp masked version of the WIYN V-band image (top) and the Ha emission (bottom). The right panels display the Ha emission and 
have been stretch to emphasize the brightest (top) and faintest (bottom) Ha emission. Grayed areas mask artifacts contaminating the Ha 
emission created by foreground stars, and some CCD artifacts. The prominent dust lane is clearly extinguishing Ha emission within the 
plane of the galaxy. When comparing the cxtraplanar filamentary-like structures in the V-band image to the Ha emission, there is little 
physical correlation between the features. The unsharp-masked Ha image, in particular, demonstrates that the few filaments found in the 
DIG are faint relative to the smooth DIG. The four faint filaments making up the large "H"-shaped structure emanating from both sides 
of the bulge region are indicated by arrows in top-right panel. The filaments originate R ~ 200 — 300 pc from the center of the galaxy and 
can be traced in our images to at least z ~ 2.2 kpc. 



can be ranked in order of decreasing DIG scale height: 
NGC 5775, NGC 891, NGC 4013, NGC 4302. This is the 
same order one would have if ranking by decreasing Ha 
intensity and decreasing prominence of DIG filaments. 

5. THE RELATIONSHIP BETWEEN 
EXTRAPLANAR IONIZED GAS AND DUST 

NGC 4013 and NGC 4302 show both dust and DIG 
in the thick disk. Is there any relationship between 
the physical structures probed by these tracers of the 
thick disk ISM? NGC 891 is the only galaxy for which 
a detailed comparison between the dust-bearing cloud 
structures and the DIG in the thick disk has been made 
at high resolution, although several works have noted 
the distinct morphologies of thick disk dust and DIG in 
NGC 891 (Dettmar 1990, HS97, HS00, Rossa & Dettmar 
2003a, 2003b, Rossa et al. 2004, THS04). Howk & Sav- 
age (HS00) argued the dusty clouds in the thick disk 
of NGC 891 trace a dense medium, in part because of 
these differing morphologies. They also estimated typ- 
ical densities for the extraplanar dust-bearing clouds of 
?ih > 2 — 10 cm -3 , significantly higher than DIG struc- 
tures at any distance from the midplane. Furthermore, 
HS00 found very few cases where structures traced in 
Ha emission and dust absorption could be even loosely 
associated. They found no direct physical connection 
between the extraplanar dust-bearing structures and the 
DIG seen in Ha emission. Rossa et al. (2004) reached the 
same conclusion after comparing the thick disk dust ex- 
tinction and Ha emission in high resolution HST imaging 
of NGC 891. Keppel et al. (1991) also noted the emission 
measure from clouds as dense as those probed by the ex- 
traplanar dust should be much higher than that observed 
from the DIG in NGC 891, suggesting the dusty clouds 
trace a different medium. 



Our observations allow us to explore the relationship 
between extraplanar dust extinction and Ha emission at 
high resolution in two additional galaxies. Extraplanar 
dust filaments are very prominent in the thick disk of 
both galaxies, whereas the extraplanar DIG layer lacks 
prominent filaments or structures in both. Figures [9] and 
[TOl show the Ha emission images on the same physical 
scales as the broadband ground based images, allowing 
for a direct comparison of the dense thick disk dust cloud 
structures and the Ha emitting DIG structures. 

In NGC 4013 we see filamentary dust structures in 
the thick disk in our broadband HST and WIYN im- 
ages along the entire radial length of the galaxy, with 
extensive and highly-structured dust filaments seen near 
the bulge region. Our images show absorbing dust struc- 
tures out to heights z ~ 2.0 kpc on both sides of the 
plane of NGC 4013. On the largest scales, the extra- 
planar dust and DIG show distinct distributions. The 
DIG is relatively smooth with weak and few filaments. 
The DIG filaments that do exist are not reflected in the 
dust. The large extraplanar dust complexes themselves 
are not directly identifiable in the Ha emission. The lack 
of a morphological connection between extraplanar dust 
clouds and the DIG is consistent between the largest and 
smallest scales in NGC 4013. 

In NGC 4302 we come to the same conclusion: the mor- 
phologies of the absorbing dust structures and Ha emit- 
ting structures show no obvious physical relationship. In 
NGC 4302 our broadband images show filamentary dust- 
bearing clouds extending to z ~ 2.5 kpc along the length 
of the galaxy. The thick disk DIG layer has very few 
filamentary structures and is smoothly distributed along 
the inner R(north) ~ 7 kpc and R(south) ~ 8.5 kpc; at 
larger R the extraplanar Ha emission quickly disappears. 
This drop-off in the DIG emission contrasts with the be- 
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Fig. 10. — LBT unsharp masked V-band image (top) and WIYN continuum-subtracted Ha image of NGC 4302 (bottom). Grayed areas 
mask artifacts due to saturated stars. The extraplanar Ha is weak in intensity except in the inner regions of the galaxy. Extraplanar 
absorbing dust structures seen in the V-band image are not reflected in the Ha emission. There is little evidence for a direct physical 
connection between the extraplanar dust, which shows structure to very small scales, and the extraplanar DIG, which is very smooth. 



havior of the extraplanar dust filaments, which are seen 
extending to large z-heights to nearly the edge of the op- 
tical disk. The Ha morphology we see in NGC 4302 is 
a result of large dust complexes extinguishing emission, 
yet no physical DIG structures are traced by any dust 
filaments seen in our HST and LBT images. Overall, 
the dust clouds are structured on very small scales and 
strongly filamentary, while the extraplanar ionized gas is 
primarily diffuse and smooth. 

When comparing the morphologies of the dust and 
DIG in the thick disk, neither NGC 4013 nor NGC 4302 
show evidence for a direct physical relationship between 
the thick disk dust seen in the broadband images and the 
extraplanar ionized gas seen in Ha emission. NGC 4302 
shows the strongest contrast in the physical appearance 
of the dust and DIG as it has one of the weakest extra- 
planar DIG layers yet has a very dusty thick disk. In a 
sense, we expect very little connection between the dust 
structures and DIG given the importance of the smooth 
components of the DIG in these two galaxies (Rand 1996, 
Rossa & Dettmar 2000, Heald et al. 2007). The smaller- 
scale structures in the dust-bearing clouds, their larger 
column densities, and even the differences in their large- 
scale distribution emphasize that the dust traces a phase 
of the thick disk ISM distinct from the DIG. 



NGC 4013 and NGC 4302 are two edge-on galaxies that 
display a lack of a direct physical correlation between 
thick disk dust and ionized gas. This physical difference 
between the morphology of the extraplanar DIG and dust 
is not unique to NGC 891. At the same time, there is a 
correlation between the presence of extraplanar dust and 
DIG disks of galaxies. As first noted by HS99 in a small 
sample of galaxies and subsequently in the much larger 
survey of Rossa & Dettmar (2003a, 2003b), when DIG 
is detected in a galaxy, extraplanar dust is also visible 
~ 90% of the time. Thus, although there is no evidence 
for a morphological relationship between the dust and 
DIG, there is a strong correlation between the presence 
of extraplanar dust and DIG in galaxies, a result reflected 
in NGC 4013 and NGC 4302. 

6. THE MULTIPHASE THICK DISKS OF SPIRAL 
GALAXIES 

We have presented deep broadband and Ha images of 
the edge-on spiral galaxies NGC 4013 and NGC 4302. 
Our images reveal significant numbers of dusty extrapla- 
nar clouds via the extinction they produce against the 
background thick disk and bulge starlight in both galax- 
ies. These galaxies were chosen in part because the mor- 
phology, brightness, and scale height of the DIG emission 
in these galaxies was significantly different than the other 
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galaxy for which a detailed comparison of the extrapla- 
nar dust and DIG is available, NGC 891 (HSOO, Rossa et 
al. 2004). Thus, our program is designed to understand 
whether the changes in the morphology and properties of 
the DIG are reflected in any way in the extraplanar dust. 
While various studies have argued that the thick disk 
dust absorption is tracing a different medium or phase 
than the DIG (Keppel 1991, HS99, HSOO, Rossa et al. 
2004) with distinct morphologies, we are testing whether 
this dense medium responds at all to the processes that 
shape the DIG properties. 

We find no evidence that the morphologies and prop- 
erties (quantity, sizes, masses, etc.) of the extraplanar 
dust-bearing clouds are subject to the strong changes like 
those seen in the DIG morphologies. Table Q] presents 
Lfir, star formation rate (SFR), and star formation 
surface density (E») for NGC 4013 and NGC 4302 
in comparison to previously studied edge-on galaxies 
NGC 5775, NGC 891, and NGC 4217. The star forma- 
tion surface density (SFR per unit area) E* = SFR/^425, 
where A25 = 7r£)| 5 /4. Ordering these galaxies by SFR 
(derived from the Lfir) Table Q] shows a connection be- 
tween the morphology of the extraplanar DIG and the 
SFR. Thus, as the importance (prominence, brightness, 
etc.) of the DIG filaments and scale height decrease 
along the sequence of NGC 5775, NGC 891, NGC 4217, 
NGC 4013, and NGC 4302, the morphological and phys- 
ical properties of the observable dust seem to vary little. 
This is in part a selection effect, since the detection of 
the extraplanar dust in our images relies on the small- 
scale structure of the absorbing structures. The thick 
disk dust is only visible because it is more opaque than 
its surroundings, and this likely implies the clouds in our 
images have significantly higher particle densities than 
their surroundings. However, this also emphasizes our 
point: the properties of the dust are disconnected from 
that of the DIG in the thick disk. 

In fact, one can see morphological and structural 
changes in the DIG that are not mirrored in the extra- 
planar dust within individual galaxies. In NGC 891 the 
southwest half of the disk has a smaller DIG scale height, 
with many fewer and less prominent filaments than the 
northeastern half of the disk (HSOO). However, there is 
relatively little change in the properties and incidence of 
dust between these regions, save that there are a very few 
DIG filaments in the northeastern half that are physically 
related to extraplanar dust structures. (Such an arrange- 
ment is only seen in this portion of NGC 891 and in some 
portions of the strongly-star forming galaxy NGC 4631; 
see Wang et al. 2001.) Similarly, the present observa- 
tions show the DIG layer in NGC 4302 cuts off sharply 
at a projected radial distance from the nucleus of R ~ 7 
- 8.5 kpc, while the extraplanar dust can be seen well be- 
yond these projected radii. Additionally, the thick disk 
near the bulge region of NGC 4013 shows highly struc- 
tured filaments whereas the extraplanar DIG morphology 
is smooth even in the brightest Ha emission. 

Apparently the factors that shape the morphology and 
quantity of the extraplanar DIG do not affect the dusty 
extraplanar clouds seen in our images. While there is a 
strong correlation between the existence of extraplanar 
DIG and extraplanar dust within a galaxy (HS99, Rossa 
& Dettmar 2003b), the ionized gas traced by Ha emis- 



sion and dense clouds traced in optical continuum ab- 
sorption have very little to do with one another. What 
drives this distinct behavior of the DIG and dense, dusty 
clouds? One possibility is that the SFR or star forma- 
tion surface density (SFR per unit area) plays a role 
in shaping the morphology of the DIG layer. Rossa & 
Dettmar (2003a) found DIG morphology was connected 
to star formation activity in a galaxy, with the transition 
from higher to lower E* or SFR being associated with a 
progression from DIG with prominent filamentary struc- 
tures through fainter, more diffuse DIG. The sequence 
then of NGC 5775, NGC 891, NGC 4217, NGC 4013, 
to NGC 4302 is one of decreasing SFR as well as de- 
creasing DIG scale height and Ha filament prominence. 
While this is a limited sample, we also see changes in star 
formation affect the DIG morphology on smaller scales. 
Rossa & Dettmar (2003a, 2003b) note that star forma- 
tion on the local and global scales were factors in the 
presence of extraplanar DIG in their survey. For ex- 
ample, the southwestern half of NGC 891 has a smaller 
DIG scale height and fewer prominent filaments than the 
northeastern half (HSOO) while also having a lower SFR 
as evidenced by radio continuum and far infrared mea- 
surements (Dahlem et al. 1994, Wainscoat et al. 1987). 
The edges of NGC 4302 beyond the end of the DIG layer 
are also representative of this trend. 

By what mechanism could the star formation or its 
surface density affect the morphology of the DIG in this 
way? The answer is unclear. It could be that the fila- 
ments seen in the DIG are tracing ionized fragments of 
supershells whose presence and brightness require very 
recent star formation in clusters to maintain the ion- 
ization, while the smoother DIG component is ionized 
by longer- lived field OB stars (e.g., Hoopes & Walter- 
bos 2000). However, Rand (1998) found that the bright 
filaments in NGC 891 did not show unusual line ratios 
(e.g., [N II]/Ha, [S II]/Ha) for the DIG. In fact, the 
filaments seemed to show the same overall behavior as 
the more diffuse DIG, with lower forbidden line to Ha 
ratios associated with brighter emission. Thus, if the 
filaments are associated with recent superbubbles, they 
still share common ionization and excitation characteris- 
tics with the general DIG. At this point this assignment 
of cause is largely speculation. It seems likely that SFR 
and/or E* likely play a significant role in not only de- 
termining whether a galaxy has significant DIG emission 
(Rossa & Dettmar 2003a, 2003b, Rand 1996), but also 
in shaping the DIG structure (Rossa & Dettmar 2003b) . 

As with the extraplanar DIG, a minimum star forma- 
tion rate or E* seems to be required before extrapla- 
nar dust is observable (Rossa & Dettmar 2003a, 2003b), 
but above that threshold there is little discernable effect 
on the morphology (Table [I}. Furthermore, there is no 
evidence that the positional changes in E* that modify 
the DIG morphology affect the dense extraplanar clouds. 
In comparing the DIG morphology sequence outlined in 
Table Q] to the dust morphologies, the DIG changes in 
connection with the SFR, while the extraplanar dust is 
both present and consistently similar between all of the 
galaxies. That the dust clouds have little associated Ha 
emission likely means they are not directly illuminated 
by ionizing radiation (Keppel et al. 1991), since they are 
optically thick to such radiation. Perhaps they have a 
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different evolutionary timescale for production than the 
DIG (e.g., Alton et al. 2000), and the stars that would 
have ionized the clouds have already died off. It may also 
be that they have a distinct radial distribution compared 
with the DIG. We see only the foreground dust-bearing 
clouds with the direct imaging techniques, but perhaps 
there are few such clouds at the smaller radial distances 
at which the DIG dominates. It is unlikely we could 
show this without velocity information for the extrapla- 
nar dust clouds, which is not likely unless they can be 
measured through CO emission (e.g., Garcfa-Burillo et 
al. 1999). 

Discussions of the origin of the matter making up the 
dense clouds seen in our images and the phase structure 
of these clouds have been presented elsewhere (HS97, 
HS99, HS00, Howk 2005). Here we expand on those dis- 
cussions briefly to reflect new thinking (both our own and 
in the literature). The interstellar thick disk clouds seen 
via their dust absorption within z ~ 2 kpc of the mid- 
planes of a number of spiral galaxies are almost certainly 
dominated by material that has been expelled from the 
thin disk via stellar feedback, perhaps both mechanical 
and radiative. In this case, the signpost used to trace 
the matter - the dust - identifies the clouds as arising 
from metal-rich material. This has two important impli- 
cations. First, the great quantities of dust in these clouds 
implies that they are dominated by metal-rich thin disk 
matter that has been expelled rather than by infalling 
(metal poor) matter as invoked to explain the change in 
DIG rotation curves with z (e.g., Fraternali & Binney 
2006, Heald et al. 2007). 

Second, the presence of dust-rich matter in the thick 
disk implies that the processes that lift the dust from 
the thin disk are not so violent that they destroy a large 
fraction of the dust. If a significant amount of dust is 
destroyed, the masses calculated for our observed clouds 
scale upwards. If only 10% of the dust survived the trip 
to the thick disk the total mass in these clouds would 
approach the total gas mass of the ISM in these galax- 
ies (HS00), which is unrealistic. The existence of dust 
in the thick disk fueled by feedback-driven expulsion is 
clearly not a new discovery. It does not cause difficulties 
with existing models for galactic fountain-type phenom- 
ena. The walls of a supershell may extend to quite large 
distances from the plane, but much of the mass was in- 
corporated into those walls while the shell expanded rel- 
atively slowly (e.g., Weaver et al. 1977) destroying little 
of the associated dust. Furthermore, it is possible that 
radiation pressure plays a role in this process (Murray et 
al. 2011, Franco et al. 1991), which would leave all of the 
dust intact. A completely independent method of prob- 
ing the existence of dust surviving far out of the plane 
is through the infrared wavelengths. NGC 891 has de- 
tectable mid-infrared radiation out to z ~ 5 kpc, which 
is dominated by emission from diffuse interstellar dust 
(Burgdorf, Ashby, & Williams 2007). Thus, extraplanar 
dust is present many kpc further from the star-forming 
disk than our images reveal. Its existence to z ~ 5 kpc 
(and even beyond; see Menard et al. 2010) implies the 
dust is not all destroyed as it is lifted from the thin disk. 

The thick disk clouds traced by the extraplanar dust 
absorption in our optical images are physically distinct 
from the DIG. They trace material with densities larger 
than that of the DIG by an order of magnitude or two. 



The large gas densities (nn > 1 — 25 cm 3 ; HS00, 
Howk 2005, this work) and high gas column densities 
(Nh > 10 20 to > 10 21 ) traced by the extraplanar dust are 
largely consistent with expectations for gas associated 
with a warm neutral or cold neutral medium. Both are 
likely present in these clouds, given the level of substruc- 
ture and evidence for core-halo type structure in some 
clouds (see also HS00). Although in principle one might 
be able to use the distribution of apparent extinction in 
isolated regions of these galaxies to disentangle the rel- 
ative contributions of these phases, the results are am- 
biguous given the fitting uncertainties in the background 
light distribution and the presence of clouds deeper in 
the galaxy with a larger fraction of unabsorbed starlight. 

HS00 and HS99 have argued that the extraplanar dust 
may to a large degree trace a thick disk CNM based on 
assessments of the physical properties and a compari- 
son of the cloud masses with Galactic molecular clouds 
(Blitz 1990, 1991). There is complementary evidence 
in our work for the presence of a thick disk CNM in 
the form of the extraplanar CO emission detected in 
NGC 4013 (Garcfa-Burillo et al. 1999) and the prob- 
able presence of thick disk H II regions tracing recent 
star formation (HS00, Rossa & Dettmar 2003a). Candi- 
date thick disk H II regions can be seen in both NGC 4013 
and NGC 4302 (point-like knots of Ha emission visible in 
Figure l9l and fTO)) . In particular, Figure [9] shows a promi- 
nent candidate to the north of the bulge of NGC 4013 
(spectroscopy of this nebula is focus of K. Rueff et al., 
in preparation). Potential areas of star formation sug- 
gest a CNM as a component of the thick disk ISM in 
these galaxies. While the masses of the clouds reported 
in Table [1] are somewhat lower than those studied in ear- 
lier works, this difference is in part due to the conscious 
choice to study smaller clouds in order to emphasize the 
small-scale structure revealed by our high resolution im- 
ages. The properties of the extraplanar dust-bearing 
structures in our images are consistent with a significant 
CNM present in the thick disks of these spiral galaxies. 
However, the fact that the clouds vastly outnumber the 
H II regions so far identified suggests that the fraction 
of the former that trace molecular clouds able to sup- 
port star formation (as opposed to "diffuse clouds" in 
the Galactic parlance) may be small. 

While the thick disk gas traced by these extraplanar 
dust structures largely originates in the thin disk, the 
production or maintenance of a thick disk CNM can come 
about in several ways. HS00 and HS99 discussed the pos- 
sibility that the CNM is consistent with standard theo- 
ries for multiphase thermal equilibrium of the ISM (e.g., 
McKee & Ostriker 1977, Wolfire et al. 1995a,b, Wolfirc 
et al. 2003). Here we add that more dynamical origins 
may be even more important. As feedback-driven out- 
flows escape from the thin disks of galaxies, they sweep 
up significant masses of material. Some of this material 
is compressed at a rate insufficient to shock the gas to 
high temperatures and/or at a time when the shell walls 
might radiatively cool. In this case, the shell serves to 
compress the gas, triggering the formation of a cold or 
even molecular medium (e.g., McCray & Kafatos 1987). 
This is, of course, the origin of the self-propagating star 
formation theories. 

In a similar vein, recent work has highlighted the in- 
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stabilities that can arise in converging flows of neutral 
material. This can serve to compress and cool the gas 
(e.g., Heitsch et al. 2008), and the convergence of the 
walls of adjacent shells may be a trigger for cloud forma- 
tion (Ntormousi et al. 2011). Our observations of extra- 
planar dust always show numerous and very complicated 
structures within the first kpc of the galaxy midplanes, 
with few readily-identifiable supershells. Some of this 
is undoubtedly due to line of sight confusion. However, 
given the difficulty of distinguishing physically connected 
structures in our images, it is not a stretch to imag- 
ine some of the complexity we see in the dust occurs 
as walls of adjacent supershells or other structures con- 
verge. Converging flow-type scenarios provide a sustain- 
able method for forming CNM clouds in the thick disk 
and for reshaping the structures in a less-ordered fashion. 
The precise mechanism of the formation of CNM clouds 
in the thick disk is little constrained by our data. How- 
ever, we note that the associations of cool (Moss et al. 
2012) or even molecular gas (Dawson et al. 2011a, b) with 
high-latitude supershells in the Milky Way offer support- 
ing evidence that the ejection of the gas from the thin 
disk may also be in part responsible for triggering the 
formation of a cool medium. The global characteristics 
and physical conditions for the material traced by the 
extraplanar dust in our images suggest the presence of 
a thick disk CNM in addition to warm neutral and ion- 
ized gas. In future works we will explore the thick disk 
star formation, perhaps from this CNM gas, implied by 
extraplanar H II regions. 

7. SUMMARY 

We have presented HST /WFPC2 and ground-based 
WIYN 3.5-m and LBT 2 x 8.4 m images of the edge-on 
spiral galaxies NGC 4013 and NGC 4302. These images 
were used to study the distribution of extraplanar dust 
and ionized gas in these galaxies. The major results of 
our work are as follows. 

1. Our new images reveal extensive distributions of 
extraplanar dust clouds in NGC 4302 and NGC 4013, 
with absorbing structures traced to heights z ~ 2 kpc 
in both galaxies. The thick disk dust-bearing gas clouds 
are seen in large complexes as well as narrow filaments, 
with significant substructure seen to the limit of our HST 
resolution (< 10 — 20 pc). 

2. We derive the properties of several extraplanar dust 
clouds in each galaxy (Table [lj. We estimate hydro- 
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gen column densities for individual dust clouds N(H) 
> 2 x 10 20 cm" 2 and masses > 1O 4 M , with the larger 
complexes having quite a bit higher masses. The implied 
particles densities range from «h > 1.0 — 4.5 cm~ 3 for 
the larger clouds, with substructure clouds having even 
larger densities, jih > 7 cm -3 . The properties of these 
individual absorbing structures lead us to conclude the 
dust clouds trace a dense phase of the thick disk ISM. 

3. Our high-resolution Ha images show the extrapla- 
nar ionized gas layers of both galaxies have a diffuse, 
smooth distribution, with very few filamentary struc- 
tures. Candidate thick disk H II regions can be seen 
in both NGC 4013 and NGC 4302. In NGC 4302, the 
Ha emission abruptly cuts off at projected radii of R ~ 7 
kpc and 8.5 kpc in the northern and southern halves of 
the galaxy, respectively. 

4. Comparing the morphologies of the extraplanar dust 
clouds to the extraplanar DIG emission we find little ev- 
idence for a direct physical connection between the dust 
and the DIG. The dust-bearing clouds are strongly fil- 
amentary while the DIG layers lack filaments and are 
smoothly distributed. The features and physical proper- 
ties of the extraplanar dust and DIG in NGC 4013 and 
NGC 4302 are distinct. 

5. The thick disk ISM extending ~ 2 kpc away from 
galactic disks is a multi-phase system, including a dense, 
cold phase traced by dust extinction and perhaps formed 
through thermal instabilities and/or in converging flows. 
This thick disk in NGC 4013 and NGC 4302 is ultimately 
fueled by galactic fountain-like expulsion from the star- 
forming thin disk that does not destroy the dust as it 
raises material to high-z. 
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TABLE 1 

General Characteristics of NGC 4013 and NGC 4302 a 



Name 


RA 


Dec. 


£>25 


V rac j 


W 20 


Dist. b 


SFR C 


Type 




[J2000] 


[J2000] 


[arcmin] 


[km s _1 ] 


[km s — 1 ] 


[Mpc] 


[Meyr- 1 ] 




NGC 4013 


11 58 31 


+43 56 48 


4.7 


831 


407 


17.0 


1.5 


Sb 


NGC 4302 


12 21 42 


+14 35 54 


4.7 


1149 


377 


16.8 


0.9 


Sc 



3 The properties presented here were taken from the NED database, except where noted. 
b We adopt distances from the Extragalactic Distance Database (Tully et al. 2009). 

c The star formation rate, SFR, were calculated following Kewley et al. (2002), using the Lpm and the Kennicutt (1998) SFR calibration 
(see SFR discussion in Sj6j . 



TABLE 2 
Log of HST/WFPC2 Observations 3 - 



Filter 


Total Exp. Time [seconds] 






NGC 4013° 


NGC 4302 


F450W (B-band) 


4000 


4000 


F555W (V-band) 


3000 


2000 


F814W (/-band) 


3000 


2000 



3 The angular and physical resolutions are O'/l and 8 pc, respectively. 

b The NGC 4013 images have non-uniform exposure times per pixel. Two separate visits with HST had distinct pointing centers and 
orientations. The central and northern sections of the galaxy, respectively, have three and two times the exposure time of the southern- 
most region. 



TABLE 3 
Log of WIYN Observations 3 



Filter 


Total Exp. Time 




Seeing b 


Resolution 




[seconds] 




[arcsec] 


[PC] 


NGC 4013 


B 


2250 




0.6 


49 


V 


2900 




0.9 


74 


I 


1800 




0.8 


66 


r 


3600 




0.8 


66 


Ha d 


14400 




0.8 


66 


NGC 4302 


r 


3600 




1.0 


82 


Ha d 


9000 




0.8 


65 


a BVI-band data were 


taken with the S2KB detector. Ha and 


r-band data were 


taken with the WTTM. 




Seeing values empirically measured in the final images. 








c Resolution based on 


assumed distances from Table [T] 








d WIYN W015 filter 














TABLE 4 








Log of LBT/LBC Observations of NGC 4302 




Filter 


Total Exp. Time 




Seeing 3 


Resolution 111 




[seconds] 




[arcsec] 


[PC] 


U 


2000 




1.2 


98 


D 


1000 




1.1 


90 


V 


1000 




1.0 


81 


I 


2000 




0.9 


73 



3 Seeing values empirically measured in the final images. 
b Resolution based on assumed distances from Table [T] 
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TABLE 5 

Properties of Individual High- z Dust Features 



Cloud a 


Dimensions 


M b 


a v c 




n H ° 


Mass f 


[J2000] 


[pcxpc] 


[PC] 


[mag.] 


[cm- 2 ] 


[cm- 3 ] 


[M ] 


NGC 4013 


115835.5+435702.0 


90 x 60 


590 


0.20 


>4xl0^ u 


2.1 


>3xl0' i 


115834.4+435713.5 


50 x 45 


640 


0.17 


^ o,,in2() 

>3x 10 


2.2 


>lxl0 4 


115833.8+435713.8 


70 x 70 


890 


0.15 


>3xl0 20 


1.3 


>2xl0 4 


115834.4+435651.9 


170 x 150 


930 


0.22 


>4xl0 20 


0.9 


>7xl0 4 


115832.9+435706.8 


40 x 30 


690 


0.22 


>4xl0 20 


4.5 


>5xl0 3 


115833.5+435648.8 


90 x 55 


850 


0.12 


>2xl0 20 


1.4 


>lxl0 4 


115833.3+435648.5 


55 x 25 


820 


0.18 


>3xl0 20 


4.4 


>lxl0 4 


115833.2+435648.4 


65 x 60 


810 


0.13 


>3xl0 20 


1.3 


>9xl0 3 


115833.1+435647.9 


110 x 70 


650 


0.21 


>4xl0 20 


1.9 


>3xl0 4 


115832.2+435644.6 


30 x 30 


690 


0.20 


>4xl0 20 


4.4 


>4xl0 3 


115832.0+435644.1 


60 x 50 


660 


0.19 


>4xl0 20 


2.3 


>lxl0 4 


115831.1+435658.2 


55 x 40 


740 


0.17 


>3xl0 20 


2.7 


>lxl0 4 


115830.3+435651.5 


45 x 40 


510 


0.18 


>4xl0 20 


2.8 


>4xl0 3 


115830.8+435640.2 


130 x 70 


500 


0.23 


>4xl0 20 


2.0 


>3xl0 4 


115832.8+435638.6 


70 x 60 


600 


0.16 


>3xl0 20 


1.8 


>lxl0 4 


115830.0+435648.9 


55 x 55 


430 


0.25 


>5xl0 20 


2.9 


>2xl0 4 


NGC 4302 


122141.9+143506.8 


90 x 80 


590 


0.14 


>3xl0" !U 


1.1 


>3xl0" 


122143.1+143512.4 


110 x 90 


900 


0.19 


>4xl0 20 


1.3 


>3xl0 4 


122143.0+143518.9 


110 x 100 


820 


0.22 


>4xl0 20 


1.4 


>6xl0 4 


122141.9+143545.4 


90 x 65 


560 


0.33 


>6xl0 20 


3.2 


>3xl0 4 


122143.1+143550.6 


60 x 50 


860 


0.26 


>5xl0 20 


3.3 


>2xl0 4 


122141.8+143601.9 


130 x 90 


660 


0.21 


>4xl0 20 


1.5 


>3xl0 4 


122143.1+143620.4 


90 x 70 


950 


0.17 


>3xl0 20 


1.5 


>4xl0 4 


122142.8+143627.3 


95 x 60 


600 


0.19 


>4xl0 20 


1.9 


>5xl0 4 


122141.8+143633.9 


95 x 65 


610 


0.31 


>6xl0 20 


2.0 


>7xl0 4 


122141.6+143649.3 


105 x 95 


820 


0.24 


>5xl0 20 


1.5 


>4xl0 4 


122143.2+143658.4 


95 x 90 


1040 


0.12 


>2xl0 20 


0.9 


>2xl0 4 



a Cloud I.D. based on R.A. and Dec. 

k Projected vertical distance from the midplanc. 

c Average apparent V-band extinction. 

Column density derived using measured ay and the Milky Way value for JVh/Av ~ 1-9 x 10 21 cm -2 mag -1 . These are lower limits since 
a v < A v . 

° Mean density derived from N H and the minor axis width. 

f Approximate mass based upon the estimated column density and projected area. Includes a factor of 1.37 correction for He. 



TABLE 6 

Comparison of Extraplanar Dust and DIG Properties 



Galaxy 




SFR b 




Dust 


DIG 


DIG 




[10 9 L ] 




[10- 3 M Q yr- 1 kpc- 2 ] 


Morphology 


Morphology 


Reference 


NGC 5775 


32.5 


9.9 


12.0 


Highly structured filaments 


Many bright filaments 


1 


NGC 891 


12.4 


3.8 


3.4 


Highly structured filaments 


Bright filaments, faint + diffuse 


2,3,4 


NGC 4217 


8.1 


2.5 


4.6 


Highly structured filaments 


Two faint patches 


5 


NGC 4013 


4.8 


1.5 


2.8 


Highly structured filaments 


Faint filaments; faint + diffuse 


6 


NGC 4302 


2.6 C 


0.9 


1.4 


Large complexes + filaments 


Faint + diffuse 


5,6 



References. — (1) Collins ct al. 2000; (2) HS00; (3) Rand, Kulkarni, & Hester (1990); (4) Dettmar (1990); (5) Rand (1996); (6) This work. 
a The Lfjr values were calculated using IRAS 60 and 100 fim fluxes from Sanders et al. (2003). 

k The star formation rates were calculated following Kewley et al. (2002) using the calculated Lfir and the Kennicutt (1998) star formation 
calibration. For NGC 5775, NGC 891, and NGC 4217 we adopt distances of 26.7 Mpc, 9.6 Mpc, and 17 Mpc, respectively. 

c The IRAS images do not resolve NGC 4302 and NGC 4298. We estimate the fractional contribution to both IRAS bands using the relative 100 
fim fluxes for these two galaxies from more recent Herschel images (Corbclli ct al. 2012). 
^ The star formation surface density, = SFR/A25, where A25 — tvD^/&. 
e Description based on images from the Hubble Legacy Archive. 



